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Abstract 

In this paper, we analytically derive an upper bound on the error in approximating the uplink 
(UL) single-cell interference by a lognormal distribution in frequency division multiple access 
(FDMA) small cell networks (SCNs). Such an upper bound is measured by the Kolmogorov-Smirnov 
(KS) distance between the actual cumulative density function (CDF) and the approximate CDF. The 
lognormal approximation is important because it allows tractable network performance analysis. Our 
results are more general than the existing works in the sense that we do not pose any requirement 
on (i) the shape and/or size of cell coverage areas, (ii) the uniformity of user equipment (UE) 
distribution, and (iii) the type of multi-path fading. Based on our results, we propose a new 
framework to directly and analytically investigate a complex network with practical deployment of 
multiple BSs placed at irregular locations, using a power lognormal approximation of the aggregate 
UL interference. The proposed network performance analysis is particularly useful for the 5th 
generation (5G) systems with general cell deployment and UE distribution. 

I. Introduction 

Small cell networks (SCNs) have been identified as one of the key enabling technologies 
in the 5th generation (5G) networks [|T]|. In this eontext, new and more powerful network per- 
formanee analysis tools are being developed to gain a deep understanding of the performanee 
implieations that SCNs bring about. These tools are significantly different from traditional 
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network performance analysis tools applicable for studying just a few macrocells only. These 
network performance analysis tools can be broadly classified into two large groups, i.e., 
macroscopic analysis and microscopic analysis [2-10]. 

The macroscopic analysis assumes that both user equipments (UEs) and base stations (BSs) 
are randomly deployed in the network, often following the homogeneous Poisson distribution, 
and usually try to derive the signal-to-interference-plus-noise ratio (SINK) distribution of 
UEs and other performance metrics such as the coverage probability and the area spectral 
efficiency [2,3]. The microscopic analysis is often conducted assuming that UEs are randomly 
placed and BSs are deterministically deployed, i.e., the BS positions are known [4-10]. 

The microscopic analysis is important because it allows for a network-specific study and 
optimization, e.g., optimizing the parameters of UL power control (81 and performing per- 
cell loading balance in a specific SCN (H . In contrast, the macroscopic analysis investigates 
network performance at a high level by averaging out all the possible BS deployments [2,3]. 
Generally speaking, the microscopic analysis gives more targeted results for specific networks 
than the macroscopic analysis, while the macroscopic analysis gives a general picture of the 
network performance. 

In this paper, we focus on the microscopic analysis. In particular, we consider an uplink 
(UE) frequency division multiple access (EDMA) SCN, which has been widely adopted in the 
4th generation (4G) networks, i.e., the UL single-carrier EDMA (SC-EDMA) system in the 
3rd Generation Partnership Project (3GPP) Long Term Evolution (LTE) networks (TTI and the 
UL orthogonal EDMA (OEDMA) system in the Worldwide Interoperability for Microwave 
Access (WiMAX) networks (T^ . Eor the UL microscopic analysis, the existing works use 

• Approach 1, which provides closed-form expressions but complicated analytical results 
for a network with a small number of interfering cells and each cell has a regularly- 
shaped coverage area, e.g., a disk or a hexagon (H. In (U, the authors considered a 
single UL interfering cell with a disk-shaped coverage area and presented closed-form 
expressions for the UL interference considering both path loss and shadow fading. 

• Approach 2, which first analyzes the UL interference and then makes an empirical 
assumption on the UL interference distribution, and on that basis derives analytical results 
for a network with multiple interfering cells, whose BSs are placed on a regularly-shaped 
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lattice, e.g., a hexagonal lattice O-Gl. Specifically, in |l5l and Q, the authors showed 
that the lognormal distribution better matches the distribution of the uplink interference 
than the conventionally assumed Gaussian distribution in a hexagonal cellular layout. 
In Ol, the authors assumed that the uplink interference in hexagonal grid based OFMDA 
cellular networks should follow a lognormal distribution. Such assumption was verified 
via simulation. 

• Approach 3, which conducts system-level simulations to directly obtain empirical results 
for a complex network with practical deployment of multiple cells, whose BSs are placed 
at irregular locations [HI, (HI, (91. In particular, the authors of dl, (HI, (9l conducted 
system-level simulations to investigate the network performance of SCNs in existing 4G 
networks and in future 5G networks. 

Obviously, Approach 1 and Approach 3 lack generality and analytical rigor, respectively. 
Regarding Approach 2, it has been a number of years since an empirical conjecture was 
extensively used in performance analysis, which stated that the UL inter-cell interference 
with disk-shaped coverage areas and uniform UE distributions could be well approximated 
by a lognormal distribution in code division multiple access (CDMA) SCNs O, @ and 
in FDMA SCNs (71. This conjecture is important since the lognormal approximation of 
interference distribution allows tractable network performance analysis. However, up to now, 
it is still unclear how accurate this lognormal approximation is. In this paper, we aim to 
answer this fundamental question, and thus making a significant contribution to constructing 
a formal tool for the microscopic analysis of network performance. Note that in our previous 
work da, we investigated an upper bound on the error of this lognormal approximation 
under the assumptions of uniform UE distribution and Rayleigh multi-path fading. In this 
paper, we will largely extend our previous work by presenting a new and tighter upper bound 
on the approximation error and remove the requirement on the types of UE distribution and 
multi-path fading. 

In this paper we focus on the analysis of UE inter-cell interference. Note that the inter¬ 
ference analysis is important because it paves way to the analyses of SINR, as well as other 
performance metrics such as the coverage probability and the area spectral efficiency. The 
contributions of this paper are as follows: 
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1) Our work analytically derives an upper bound on the error in approximating the UL 
single-eell interference in FDMA SCNs by a lognormal distribution. Sueh error is mea¬ 
sured by the Kolmogorov-Smimov (KS) distance [[T3ll between the aetual eumulative 
density function (CDF) and the approximate CDF. 

2) Unlike the existing works on the mieroseopie analysis, e.g., [4-11], our work does 
not pose any requirement on (i) the shape and/or size of eell eoverage areas, (ii) the 
uniformity of UE distribution, and (iii) the type of multi-path fading. Thus, our proposed 
framework is more general and useful for network performanee analysis. 

3) Based on our work, a new approach can be established to fill an important theoretical 
gap in the existing mieroseopie analysis of network performance, which either assumes 
very simple BS deployments or relies on empirical results. Sueh new approaeh allows 
us to directly investigate a complex network with practical deployment of multiple 
BSs placed at irregular locations, while retaining mathematical rigor in the analysis. 
In order to do that, we first verify the aeeuracy of the approximated UL interference 
distribution for eaeh small eell, then we approximate the aggregate UL interferenee by 
a power lognormal distribution. Speeifically, the CDL of a power lognormal distribution 
is a power funetion of the CDL of a lognormal distribution. 

The remainder of the paper is structured as follows. In Seetion HB the network seenario 
and the system model are described. In Section |nll our approach to studying the UL inter¬ 
cell interferenee in LDMA SCNs is presented, followed by the validation of our results via 
simulations in Seetion |IV[ Linally, the eonclusions are drawn in Seetion |V1 

II. Network Scenario and System Model 

In this paper, we consider UL transmissions and assume that in one frequeney resource 
bloek (RB) and in a given time slot, only one UE is scheduled by eaeh small cell BS to 
perform an UL transmission, whieh is a reasonable assumption in line with the 4G networks, 
i.e., the UL SC-LDMA system in the LTE networks ifTTI and the UL OLDMA system in 
the WiMAX networks [fT^ . We assume that eaeh small cell has at least one assoeiated UE 
beeause the small eell BSs having no UE do not contribute to the uplink interference analyzed 
in this paper, thereby can be ignored. 
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Regarding the network seenario, we eonsider a SCN with multiple small eells operating 
on the same earner frequeney as shown in Fig. [TJ In Fig. [H the SCN eonsists of B small 
cells, each of which is managed by a BS. The network includes one small cell of interest 
denoted by Ci and B — 1 interfering small cells denoted by Ci,,b G { 2 , ..., B}. We focus 
on a particular frequency RB, and denote by K^, the active UE associated with small cell Cb 
using that frequency RB. Moreover, we denote by Rb the coverage area of small cell Cb, in 
which its associated UEs are randomly distributed. Note that the coverage areas of adjacent 
small cells may overlap due to the arbitrary shape and size of { Rb } ,b G [ 2 , ..., B}. 



Fig. 1. A schematic model of the considered SCN. 

The distance (in km) from the BS of Cb to the BS of Ci, b E {1,..., 5}, and the distance 
from UE Kb to the BS of C^, b,m E {1,..., S}, are denoted by Db and dbm, respectively. In 
this paper, we consider a deterministic deployment of BSs, and thus the set {Db} is assumed 
to be known. However, UE Kb is assumed to be randomly distributed in Rb with a distribution 
function (z), z E Rb- Hence, dbm is a random variable (RV), whose distribution cannot 
be readily expressed in an analytical form due to the arbitrary shape and size of Rb, and the 
arbitrary form of fz^ (z). Unlike the existing works, e.g., [2-8], that only treat uniform UE 
distributions, in this work we investigate a general probability density function (PDE) of UE 
distribution denoted by fz^ (z), which satisfies 0 < fz^ (z) < +oo, 2 G Rb, and its integral 
over Rb equals to one, i.e., fz^ {z) dz = 1. 

In the following, we present the modeling of path loss, shadow fading, UE transmission 
power, and multi-path fading, respectively. 

Based on the definition of dbm, the path loss (in dB) from UE Kb to the BS of Cm is 
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modeled as 

Lbm ^ “ 1 “ ( 1 ) 

where A is the path loss at the reference distance of dbm = 1 and a is the path loss exponent. 
In practice, A and a are constants obtainable from field tests lfT4ll . Note that Lbm is a RV 
due to the randomness of dbm- 

The shadow fading (in dB) from UE Kb to the BS of Cm is denoted by Sbm, and it is usually 
modeled as a zero-mean Gaussian RV because the linear-scale value of Sbm is commonly 
assumed to follow a lognormal distribution [[T4l . Hence, in this paper, we model Sbm as an 
independently and identically distributed (i.i.d.) zero-mean Gaussian RV with variance 
denoted by Sbm ~ AC (0, o-|had) ■ 

The UL transmission power (in dBm) of UE Kb is denoted by Pb. In practice, Pb is 
usually subject to a semi-static power control (PC) mechanismQ, e.g., the fractional pathloss 
compensation (EPC) scheme lfT4]l . Based on this PPG scheme, Pb in dBm is modeled as 

Pb = Po A V (Lbb + Sbb) , (2) 

where Pq is the target received power at the BS in dBm on the considered frequency RB, 
r] G (0,1] is the EPC factor, Lbb has been defined in dU), and Sbb ~ AA(O)^shad) has been 
discussed above. 

The multi-path fading channel from UE Kb to the BS of Cm is denoted by hbm £ C, where 
we assume that each UE and each BS are equipped with one omni-directional antenna. In this 
paper, we consider a general type of multi-path fading by assuming that the effective channel 
gain (in dB) associated with hbm is defined as lOlogj^g l^bml and denoted by Hbm, which 
follows an i.i.d. distribution with the PDE of fn {h). Eor example, \hbm\^ can be respectively 
characterized by an exponential distribution or a non-central chi-squared distribution in case 
of Rayleigh fading or Rician fading [fTSlI . 

'Note that in practice Ph is also constrained by the maximum value of the UL power, denoted by Pmax at the UE. 
However, the power constraint is a minor issue for UEs in SCNs since they are generally not power-limited due to the 
close proximity of a UE and its associated SCN BS. For example, it is recommended in 1141 that Pmax is smaller than the 
SCN BS downlink (DL) power by only IdB, which grants a similar outreach range of signal transmission for the BS and 
the UE. Therefore, the UL power limitation is a minor issue as long as the UE is able to connect with the BS in the DL. 
For the sake of tractability, in this paper, we model Ph as ill, which has been widely adopted in the literature [3,6-8,11]. 
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III. Analysis of the UL Interference Distribution 

Based on the definition of RVs diseussed in Seetionllll the UL received interference power 
(in dBm) from UE Ki, to the BS of Ci ean be written as 

- Lw - + i/w 


— -Pq + ivLfjb — Lbi) + {rjSbb — Pfei) + Hbi 

= {Pq + L + S') + Hbii 

= I^^+Hb^, (3) 


where dll) is plugged into the step (a) of dH), and L and S are defined as L = [rjLbb — Lbi) 
and S = {rjSbb — Sbi), respeetively. Apparently, L and S are independent RVs. Besides, the 
first part of h is further defined as = (Pq + L + S). Sinee Sbh and Sbi (6 G {2,..., B}) 
are i.i.d. zero-mean Gaussian RVs, it is easy to show that S is also a Gaussian RV, whose 
mean and varianee are 


{ds = 0 

|cr| = (1 + 7?^) alhad 


(4) 


From the definition of h in dS), the aggregate interference power (in mW) from all 
interfering UEs to the BS of Ci ean be written as 

B 

= (5) 

6=2 

In the following subseetions, we will analyze the distribution of in three steps: 

• First, we investigate the distribution of shown in dU) and approximate sueh distribu¬ 
tion by a Gaussian distribution. The upper-bound of the approximation error measured 
by the Kolmogorov-Smimov (KS) [[T3ll distanee is derived in elosed-form expressions. 

• Seeond, we analyze the sum of RVs -f Hbi^ shown in d3]) and further approximate 
the distribution of h by another Gaussian distribution. The upper-bound of the approx¬ 
imation error measured by the KS distanee lfT3]l is derived in elosed-form expressions. 

• Third, we show that the distribution of ean be well approximated by a power 
lognormal distribution, i.e., with the CDF being a power funetion of the CDF of a 
lognormal distribution. 
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A. The Distribution of in (13) 

Considering the complicated mathematical form of L in ([3]), we can find that the PDF of 
L is generally not tractable because L is a RV with respect to d^b and dbi, which jointly 
depend on the arbitrary shape/size of Rb and the arbitrary form of fzi, (z). Specifically, for 
any point z e Rb, dbb and dbi are geometric functions of 2 ;, and the probability density of 2 : is 
(z). Despite the intractable nature of L, we will show that can still be approximated 
by a Gaussian RV with bounded approximation errors. In more detail, we investigate an 
upper bound on the error in approximating the sum of an arbitrary RV and a Gaussian RV, 
i.e., (L + S), by another Gaussian RV. To that end, we denote by and the mean and 
variance of L respectively. Moreover, we define two zero-mean RVs L = L — and 
S = S — respectively. As a result, in (|3]) can be re-written as 

dj, ^ ~ (^L + s'j + {Pq + fiL + bs) ■ ( 6 ) 

Next, we approximate (^L + by a Gaussian RV. And it follows that can also be 
approximated by the same Gaussian RV with an offset (Pq + hi + Ts)- 

1) The Distribution of {ju + in 

For the convenience of mathematical expression, the mean, the variance, the 3rd moment 
and the 4th moment of L are denoted by = 0, cr| = cr|, and , respectively. 
Besides, considering dH), the mean and the variance of S are denoted by = 0, (t| = (t|, 
respectively. Moreover, we define G = L + S and denote the PDF of L by fp (/). 

In order to quantify the approximation error between the distribution G and its approximate 
Gaussian distribution, we invoke the following definition on the Kolmogorov-Smimov (KS) 
distance between two CDFs ifTSll . 

Definition 1. Suppose that the CDFs of RVs X and X are Fx {x) and (x), respectively. 
Then the Kolmogorov-Smimov (KS) distance between Fx (x) and (x) is defined as 

KS (X,Y) = sup jFx (x) - (x)j. (7) 

The KS distance is a widely used metric to measure the difference between two CDFs. 
Based on Definition [H we present Theorem [21 in the following to bound the KS distance 
between the CDF of G and that of the corresponding approximate zero-mean Gaussian RV 
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with a variance of (cri, + cr|). 


Theorem 2. Considering the zero-mean RV G = L-\-S given by (HI), the KS distance between 
the CDF of G and that of the corresponding approximate zero-mean Gaussian RV with a 
variance of {a\ + a‘g) is bounded by 


sup 




V^L + 


(Tt 


< max{ei + 82 , £ 3 } , 


( 8 ) 


where Fq (g) and $ (■) are respectively the CDF of G and the CDF of the standard normal 
distribution, and max {x,y} extracts the larger value between x and y. 


Furthermore, Ei in (jS]) is expressed as 

1. . \ ^ X ( (^1 + h ) + 


y/2as 


(Tq \ 1 ^ I u 

' + 


ki 


C^s 


,.,^1,(9) 


where u, p, ki and k 2 are positive scalars. Besides, fi (/) is the PDF of L, and 5 q {u,p,k) 
and (5i {u,p, k) are given by 

5o{uj,p, k) = -;^erfc{{ 2 p+l)uj) + erfc{f^ - k) 

,<^1 = 7l:^''/c((2p+ 1)0;) + 2 

where erfc{-) is the complementary error function HTSll . 

Moreover, 82 in (jS]) is given by 


( 10 ) 


£2 = 


TT 


2p-l 

n=l,n odd 


where Vn = - exp 


2 iPl[-Ts)’ ^n = iexp 


rPu)'^ ( 

I 


( 11 ) 

and Pi {t) is the 


characteristic function / I7^ of fi (/). 
Finall, 83 in ([S]) is expressed as 


£3 = ^ + ]^erfc{ki] 


( 12 ) 


Proof: See Appendix A. ■ 

Theorem [2] is useful to quantify the maximum error of approximating G by a Gaussian 
RV. From the proof in Appendix A, it ean be seen that 

• £1 in dH]) is eaused by the residual errors from the p-truncated Fourier series expansion 
of erfe (■), where p is the number of the truneated terms. 
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• £2 in dH) is the major contributor to the derived upper-bound of the KS distance in 
Theorem [2l 

\g\ > /ciwhere ki takes a large value to show the asymptotic behavior of the 
interested CDFs. 

To obtain some insights on the typical values of Si, £2 and £3, we should first discuss the 
appropriate choices of the values of u, p, ki and k2. To that end, we have the following 
remarks. 

Remark 1: ce is the fundamental frequency of the Fourier series expansion of erfc (■), and 
thus it should satisfy the condition that ^ is much larger than the spread of x in erfc (x) [[Bl, 
which can be estimated to be in the range of [—5,5] because |erfc (—cxo) — erfc(—5)| = 
1.54 X 10“^^ ~ 0 and |erfc (+cxo) — erfc(5)| = 1.54 x 10“^^ ~ 0. Therefore, ^ should be 
much larger than 10 due to the range of [—5, 5]. We propose that a typical value of ^ could 
be 1000, i.e., oj = 0 .001. 



£3 in dH) measures the asymptotic difference between Fg (g) and $ 


Remark 2: p is the number of the truncated terms in the p-truncated Fourier series 
expansion of erfc (•), and thus it should be sufficiently large to make the residual error 
caused by the truncation of the Fourier series expansion very small. The part of the residual 
error that is related to p, is expressed as 5 {uj,p) = ;^^erfc ((2p + 1) u) in (fT^ . From this 
expression, we can see that p should be larger than ^ to make S{uj,p) sufficiently small, 
because erfc ((2 x ^ + l) w) = 1.4 x 10“® when u = 0.001. We propose that a typical value 
of p could be which equals to 4000 when u = 0.001. The corresponding 6 {u,p) will then 
drop to a very small value of 1.1 x 10“^^. 

Remark 3: Since ki is related to the asymptotic behavior of Fg {g) and $ ( . ^ 

W^L+^S 

when |g'| > ki^^a\ + (t|, and thus ki should be sufficiently large to make £3 a small value. 
We propose that the typical value of ki should be at least 100 to achieve an £3 as small as 
1.0 X 10“^, which corresponds to a small error of ±0.01 percentile. As will be shown in 
Section |I3 ki is set to 500, and we have £3 = 4 x 10“®, which is confirmed to be always 
smaller than £2 in our numerical results. 

Remark 4: As discussed in Appendix A, the introduction of /c 2 is to facilitate the bounding 
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of the integral over I from —cxo to +cxo. From Q, the first term of (HI), i.e., \5i {uj,p) 
dominates ei because Jq k) is a very small value with appropriate choices of co, p, ki 
and /c 2 . For example, as discussed above, we choose u = 0.001, p = 4000, ki = k 2 = 500, 
and CTx, = as, then the first term of (HI) becomes 4 x 10“®, while the sum of the rest terms 
is merely 6.3 x 10“^^. Since a larger k 2 directly leads to a smaller Si, we propose that the 
typical value of /c 2 should be at least 100 to achieve an ei as small as 1.0 x 10 “^, which 
corresponds to a small error of ±0.01 percentile. As will be shown in Section |I3 k 2 is set 
to 500, and we have ei = 4 x 10“®, which is confirmed to be always smaller than 82 in our 
numerical results. 

Considering the expressions of and £3 shown in (HI) and (fT^ . respectively, we investigate 
a special case of ki = k 2 and it is straightforward to get £1 > £ 3 . Hence, we propose the 
following corollary to simplify Theorem [2l 

Corollary 3. When ki = k 2 , the KS distance between the CDF of G and that of the 
corresponding approximate zero-mean Gaussian RV with a variance of {a\ ± cr|) is bounded 


by 



(13) 


Note that compared with ([ 8 ]), £3 does not exist in the right-hand side of (fTSl) . which largely 
simplifies our analysis and discussion. Therefore, in the sequel we will only consider (fT3]) to 
quantify the maximum error of approximating G* by a Gaussian RV. 

With the discussed typical values of u, p, ki and /c 2 , we can see that £1 in (fT3l) can 
be controlled to be as small as 4 x 10“®, which leaves £2 as the major contributor to the 
derived upper-bound of the KS distance. We will briefly discuss the calculation of £2 in the 
next subsection. Note that further optimization of ki and k 2 to reduce £3 and £1 even below 
4 X 10“® is possible. However, such optimization has a marginal impact on the derived upper 
bound of the approximation error because £2 is independent of ki and /c 2 . 

With Theorem [2] and Corollary [3] characterizing the upper bound of the approximation 
error, we propose to approximate in (© by a Gaussian RV G^, whose mean and variance 
can be computed by 
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{ fJ'Gb — Po + fJ-L + fl'S 

^Gb =(^L + (^S 


(14) 


2) The Calculation of 82 in m-- 

For each small cell Cb, considering the definition of L and L respectively presented in ([3]) 
and ®, we can evaluate TL small cell Cb as 


TL 


■—)= r°°exp (-j—)h{l)dl 

/ J-00 V 0-5 / 

f f ■ L{z) - fiL\ ^ f 

= exp -jnu - fzb (z) dz 

JRb \ (^S J 

(a) f f . T]Lbb{z) - Lbi{z) - hl\ r / X , 

= J exp - — - ) fzb [Z) dz, 


(15) 


where Q should be inserted into the step (a) of (fTSl) and /i^ is the mean of L. 

With the result of ^, we can then compute 82 according to its definition in (fTTI) . 


B. The Distribution of Ib in (13) 

Having approximated by a Gaussian RV Gb, we can approximate ([3]) as 

Ib ~ Gb + Hbi- (16) 

It is interesting to note that, similar to in ([31), the approximate expression of h in (fT^ 
also contains a Gaussian RV Gb and an arbitrary RV Hbi with the PDF of fn (h). 

Based on the above observation, we propose to reuse Theorem[2]and Corollary [3] to quantify 
the error in approximating f in (fT^ by an another Gaussian RV. To that end, similar to I©, 
we define two zero-mean RVs as Gb = Gb — fJ-Gb ^nd Hbi = Hbi — hHbi ^ where ^Cb ^Hbi 
are the means of Gb and Hbi, respectively. Besides, the variance of Gb and Hbi are denoted 
by and , respectively. As a result, we can re-formulate (fT^ as 

h ~ + Hbl^ + ihGb + hHbi) ■ (17) 

Next, we approximate + Hbi^ by a Gaussian RV. And it follows that Ib can also be 
approximated by the same Gaussian RV with an offset (^^^ + fCHbi)- 

1) The Distribution of (^b + in i lTZl) .- 
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We propose to approximate + Hbij as a zero-mean Gaussian RV with a variance 
of aQ^ + , then from Corollary [3l the error measured by the KS distance between the 

actual CDF and the approximate CDF can be upper-bounded by e\ + e'2, where e\ and e'2 are 
respectively computed using (|9l) and (fTTI) with the following RV changes, 

(18) 




S —)■ Gb, S —)■ Gb, L —)■ Hbi, L —)■ Hbi 

[4 = 

£2 

S —)■ Gb, S —)■ Gb, L —)■ Hbi, L — 7 - Hbi 


where X —> V denotes the RV change of replacing X with Y. 

2) The approximate PDF and CDF of h: 

With the approximation of {Ch + Hbi^ as a zero-mean Gaussian RV, we propose to 
approximate 4 in (fT71) as another Gaussian RV Qb, whose mean and variance are 

dQb = TCb + dHbi 


(19) 




According to Theorem [2] and Corollary [3l and from ® , (fTTI) and (fTSf), the total error of 
approximating Ib as Qb, measured by the KS distance between the CDF of h and that of 
Qb, can be upper-bounded by e as 

e = [si + 82) + (^1 + 82) , ( 20 ) 


where (ei -I- 82 ) is the upper bound on the error associated with the approximation of 
in (© as a Gaussian RV Gb, and { 8 \ + 82 ) is the upper bound on the error associated with 
the approximation of Gb + Hbi in (IT^ as a Gaussian RV Qb- 


C. The Distribution of in (fj]) 

The study on the approximate distribution of the sum of multiple independent lognormal 
RVs has been going on for more than five decades IfTTfl - ETI . According to IfTTlI - lfT^ . the 
sum of multiple independent lognormal RVs can be well approximated by another lognormal 
RV However, some recent studies lfT9ll - [l^ concluded that the sum of multiple independent 
lognormal RVs is better approximated by a power lognormal RV, i.e., with the CDF being 
a power function of $ (■)• In this paper, we adopt the power lognormal approximation for 
jniw, which will be explained in the following. 
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In our case, since each Ib,b G {2,..., B} is approximated by the Gaussian RV Qb, the 
sum of can be well approximated by a power lognormal RV [[T^ - [|2T]| expressed as 

jmw _ where the PDF and CDF of Q can be respectively written as IfT^ 


PDF of g : /q (g) = A$ 
CDF of g : Fq (g) = $ 






^exp 
Q 




( 21 ) 


- rhA I 

o'Q 


where the parameters A, /rg and ctq are obtained from and {cq^}- The method to 

accomplish such task has been well addressed in [[T9ll - ll2T]l . In Appendix B, we provide an 
example to obtain A, fig and ag based on IfTTlI . If20]| and ll2TII . 

As a result of (ED, the PDF and CDF of can be respectively written as 
jpDFof/-: = ^ 

I CDF of (t;) = 

where ( = is a scalar caused by the variable change from 10 log^ v to In v. 

Finally, we propose that the distribution of J™"" can be approximated by that of shown 
in (l22l) . Note that in this step of approximation, the approximation error is dependent on the 
adopted approximate distribution of the sum of multiple independent lognormal RVs. We will 
study such approximation error in our future work. Note that some recent studies IfT^ - ETI 
have shown that the error associated with the power lognormal approximation is reasonably 
small and good enough for practical use. 


D. Summary of the Proposed Analysis of the UL Interference Distribution 

To sum up, in the following, we highlight the main steps in our proposed microscopic 
analysis of the UL interference distribution. First, for each b G {2,..., B}, we use (0 
and (HTI) to check (^i + £ 2 ) associated with the approximation of in ® as a Gaussian 
RV Gb- Second, for each b G {2,..., B}, we use (fTSl) to check + e' 2 ) associated with 
the approximation of Gb + Hbi in (fT^ as a Gaussian RV Qb- The upper bound of the 
total approximation error of the above two steps is further obtained from (EOl) . without any 
requirement on (i) the shape and/or size of cell coverage areas, (ii) the uniformity of UE 
distribution, and (iii) the type of multi-path fading. Finally, we approximate in dD by a 
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power lognormal RV 1'^'^ = lOio'?, where the PDF and the CDF of Q is expressed as (l2TI) 
with the parameters A, /ig and ag obtained from, e.g., Appendix B. 

IV. Simulation and Discussion 

In order to validate the approximation from the proposed mieroscopic analysis of the UL 
interferenee, we eonduet simulations considering two types of scenarios, i.e., one with a single 
interfering cell and the other with multiple interfering cells. As discussed in Remarks 1~4 
in Subsection IIII-A1 1 the parameters for evaluating Theorem and Corollary [3] are set to: 
u = 0.001, p = 4000, ki = k 2 = 500. According to the 3GPP standards HU, the system 
parameters are set to: A = 103.8, a = 20.9, Pq = —76 dBm, p = 0.8, and as = 10 dB. 
Besides, the minimum BS-to-UE distance is assumed to be 0.005 km [[T4ll . 

A. The Scenario with a Single Interfering Cell 

In this scenario, the number of BSs B is set to 2, and UEs in small cell C 2 are assumed to 
be uniformly distributed and non-uniformly distributed in the coverage area R 2 , respectively. 
1) Uniformly Distributed UEs: 



In this case, we consider uniformly distributed UEs, as shown in Eig. [21 The x-markers 
indicate BS locations where the BS location of Ci has been explicitly pointed out. The dash- 
dot line indicates a reference disk to illustrate the reference size of small cell C 2 . The radius 
of such a reference circle is denoted by r, and the distance between the BS of Ci and the BS 
of C 2 , i.e., D 2 , is assumed to be 1.5r. Note that 1.5r is just an example of D 2 and the specific 
value of D 2 has no impact on the procedure of our analysis. In our simulations, the values 
of r (in km) are set to 0.01, 0.02 and 0.04, respectively [HI. In this scenario, the interfering 
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UE K 2 is uniformly distributed in an irregularly shaped coverage area R 2 , as shown by the 
area outlined by the solid line in Fig. [2l The shape of R 2 is the intersection of a square, 
a circle and an ellipse, which has a very complicated generation function. Examples of the 
possible positions of K 2 within R 2 are shown as dots in Fig. 

Despite the complicated shape of R 2 , our proposed microscopic analysis of the UE inter¬ 
ference distribution still can be applied. Specifically, from Corollary [3l dH), (fTTl) . and (fTSl) . we 
can examine the validity of approximating R as a Gaussian RV by checking the corresponding 
e given by (l20l) . The results of e are tabulated for various values of r with the assumption 
of Rayleigh fading in Table HI Moreover, e is broken into (ei -f 62 ) and {e\ + € 2 ), with the 
marginal contribution of ei and shown in Table U as well. From this table, we can observe 
that 82 and e '2 are about 100 times larger than ei and ef indicating the dominance of 82 and 
e '2 in the result of £. For all the investigated values of r, the values of e are below 0.01. 
Consequently, the approximation of R as a Gaussian RV should be tight. Note that r = 0.01 
and r = 0.02 correspond to the typical network configurations for future dense and ultra- 
dense SCNs [[n, which shows that our proposed microscopic analysis of the UE interference 
distribution can be readily used to study future dense and ultra-dense SCNs. Also note that 
the upper bound e is much tighter than that presented in our previous work ifTOll . Specifically, 
when r = 0.04, the upper bound on the approximation error in [fTOll is around 1.84 x 10“^, 
while in Table |I]e is only 4.6 x 10“^, which is a significant improvement compared with our 
previous result in ifTOll . 


Table I 

Approximation errors of the proposed analysis 
{B = 2, UNIFORMLY DISTRIBUTED UES, RAYLEIGH FADING). 


r 

£ 

El + £2 

£1 

e'l + £2 

e'l 

MQ2 

6^02 

Actual error 

0.01 

4.9 X 10-3 

3.5 X 10-4 

4 X 10-6 

4.6 X 10-3 

4 X 10-6 

-97.1 

205.3 

1.9 X 10-3 

0.02 

5.0 X 10-3 

5.2 X 10-4 

4 X 10-6 

4.5 X 10-3 

4 X 10-6 

-99.7 

207.7 

1.9 X 10-3 

0.04 

4.6 X 10-3 

1.1 X 10-4 

4 X 10-6 

4.5 X 10-3 

4 X 10-6 

-101.5 

211.4 

1.8 X 10-3 


To further verify the accuracy of the proposed approximation, we plot the simulation results 
of the CDF of R and the analytical results of the approximate Gaussian CDF according to 
(fT9l) for the considered R 2 in Fig. [3l The numerical results of and aq^ are also listed 
in Table HI As can be seen from Fig. [3l the proposed Gaussian approximation of R is very 
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tight for the considered i ?2 with such an irregular shape. Note that the actual errors between 
the simulation and the approximation of the CDF of I 2 are also shown in Table IH the results 
of which establish the validity of the derived upper bound on the approximation error. 



Fig. 3. The simulation and the approximation of the CDF of I 2 (B — 2, uniformly distributed UEs, Rayleigh fading). 

As discussed above, Table H] is obtained with the assumption of Rayleigh fading. In the 
following, we will check the approximation errors for Rician fading and another case without 
multi-path fading. According to IfTSlI . a ratio between the power in the line-of-sight (LOS) 
path and the power in the other scattered paths should be defined for Rician fading. Such 
ratio is denoted by T in our paper. Note that Rician fading will degrade to Rayleigh fading 
when r = 0. For various values of r under the assumption of Rician fading with F = 10, 
the results of e are tabulated in Table HB Comparing Table HJ with Table H we can see that 
the difference lies in the values of {e'l + e' 2 ) because it solely depends on the assumption 
of multi-path fading as discussed in Subsection ITTT-R1 1 Note that the errors caused by the 
consideration of Rician fading is actually smaller than those of Rayleigh fading, because 
Rician fading incurs less randomness than Rayleigh fading due to the dominant LOS path 
component in the multi-path fading. Such reduction of randomness can be further observed 
in an extreme case of F = -|-cxo, where a deterministic LOS path completely dominates the 
multi-path fading. The results of e for such extreme case with F = -|-cxo, i.e., no multi-path 
fading, are exhibited in Table Hill where and are all zeros because the approximation 
step addressed in Subsection IIII-Bl is skipped, rendering Hi,i = OdB. 

For all the investigated values of r, the values of e are smaller than 0.001 in both Table HI] 
and Table nm indicating that the approximation of I 2 as a Gaussian RV should be even tighter 
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than that observed in Fig. [3l For brevity, we omit these figures. Note that the actual errors 
between the simulation and the approximation of the CDF of I 2 are provided in Table HI] and 
Table Hill the results of which also establish the validity of the derived upper bound e. 


Table II 

Approximation errors of the proposed analysis 
{B = 2, uniformly distributed UEs, Rician fading with T = 10). 


r 

£ 

Si + £2 

£1 

e[ -1- £2 

£'i 

MQ2 


Actual error 

0.01 

5.4 X 10-^ 

3.5 X 10-4 

4 X 10-® 

1.9 X 10-4 

4 X 10-6 

-95.0 

178.2 

2.0 X 10-4 

0.02 

7.0 X 10-4 

5.2 X 10-4 

4 X 10-6 

1.9 X 10-4 

4 X 10-6 

-97.6 

180.8 

4.8 X 10-4 

0.04 

3.0 X 10-4 

1.1 X 10-4 

4 X 10-6 

1.9 X 10-4 

4 X 10-6 

-99.4 

184.5 

1.9 X 10-4 


Table III 

Approximation errors of the proposed analysis 
{B = 2, UNIFORMLY DISTRIBUTED UEs, NO MULTI-PATH FADING). 


r 

£ 

£1 + £2 

£1 

£1 + £2 

£'i 

MQ2 

to 

Actual error 

0.01 

3.5 X 10-4 

3.5 X 10-4 

4 X 10-6 

0 

0 

-94.6 

174.2 

1.3 X 10-4 

0.02 

5.2 X 10-4 

5.2 X 10-4 

4 X 10-6 

0 

0 

-97.2 

176.7 

2.4 X 10-4 

0.04 

1.1 X 10-4 

1.1 X 10-4 

4 X 10-6 

0 

0 

-99.0 

180.5 

0.6 X 10-4 


2) Non-Uniformly Distributed UEs: 

In this subsection, we further investigate the scenario discussed in Subsection IIV-A1[ 
Different from previous assumption, here we consider that the interfering UE K 2 is no 
longer uniformly distributed in R 2 . In this subsection, we consider a UE distribution function 
expressed as fz2 {z) = ^,z e R2, where p is the radial coordinate of z in the polar 
coordinate system, the origin of which is placed at the position of the BS of C 2 . Besides, W 
is a normalization constant to make fz 2 (z) dz = 1. In the considered non-uniform UE 
distribution, UEs are more likely to locate in the close vicinity of the BS of C 2 , as shown in 
Eig. m where examples of the possible positions of K 2 within R 2 are shown as dots. Note 
that the considered fz 2 (z) is just an example of the non-uniformly distributed UEs in R 2 , 
which reflects a reasonable network planning that BSs have been deployed at the center spots 
of UE clusters. Other forms of the UE distribution function will not affect the procedure of 
our analysis, only the approximation error values may change with the choice of the UE 
distribution function. Since we have shown in Subsection IIV-AII that Rayleigh fading is the 
worst case for the proposed analysis, we adopt the assumption of Rayleigh fading in this 
subsection. 
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From (l20l) . we can evaluate the quality of approximating I 2 as a Gaussian RV by checking 
the corresponding e. Like Table H the results of e for this network scenario are tabulated for 
various values of r in Table |I3 From Table |I3 we can see that the values of e are small, 
i.e., below 0.01, which indicates that the approximation of I 2 as a Gaussian RV should be 
tight, as can be confirmed from the actual error values in Table |IV1 


Table IV 

Approximation errors of the proposed analysis {B = 2, non-uniformly distributed UEs with the 
DISTRIBUTION FUNCTION OF fz^ (z) = ^,Z £ R 2 , RAYLEIGH FADING). 


r 

£ 

£1 + £2 

£1 

£1 + £2 

A 


6^02 

Actual error 

0.01 

4.9 X 10-3 

4.0 X 10-4 

4 X 10-6 

4.5 X 10-3 

4 X 10-6 

-97.13 

205.1 

1.8 X 10-3 

0.02 

5.3 X 10-3 

7.1 X 10-4 

4 X 10-6 

4.6 X 10-3 

4 X 10-6 

-100.6 

207.9 

2.2 X 10-3 

0.04 

4.9 X 10-3 

4.0 X 10-4 

4 X 10-6 

4.5 X 10-3 

4 X 10-6 

-103.2 

214.4 

1.7 X 10-3 


B. The Scenario with Multiple Interfering Cells 

In this subsection, we apply the proposed framework to a more complex network with 
practical deployment of multiple cells and provide the approximation of the UL interference 
distribution. 



Fig. 5. Illustration of a practical deployment of multiple small cells (r = 0.04 km). 
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Here, we eonsider a 3GPP-eompliant scenario lfT4l . as shown in Fig. [51 where the number 
of BSs B is set to 84 and all small cell BSs are represented by x-markers. Particularly, the 
BS of Cl has been explicitly pointed out. The reference coverage area for each small cell is 
a disk with a radius of r mi. As in previous subsections, the values of r (in km) are set 
to 0.01, 0.02 and 0.04, respectively. The reference disk-shaped areas can be easily seen in 
Fig. [5] from any isolated small cell. However, due to the irregular positions of the cells, the 
actual coverage areas of the considered cells are of irregular shapes. The irregularly shaped 
coverage areas are outlined by solid lines in Fig. [51 

An important note is that the considered network scenario is different from that adopted 
in [2,3], where coverage areas are defined as Voronoi cells generated by the Poisson dis¬ 
tributed BSs and the entire network area is covered by those Voronoi cells. In practice, small 
cells are mainly used for capacity boosting in specific populated areas, rather than the provi¬ 
sion of an umbrella coverage for all UEs [[U. In this light, the 3 GPP standards recommend 
the hotspot scenario depicted in Fig. [5] for UE distribution in the performance evaluation 
of practical SCNs and we adopt such network scenario in this subsection. Nevertheless, we 
should mention that the proposed microscopic analysis of the UE interference distribution 
can still be applied on a particular Voronoi tessellation. This is because Theorem [2] and 
Corollary [3] in this paper do not rely on particular shape and/or size of coverage areas. 

In this subsection, we consider the same UE distribution function as the one discussed in 
Subsection IIV-A2[ i.e., fz^ {z) = G Rf,, where p is the radial coordinate of z in the 
polar coordinate system with its origin placed at the position of the BS of Ch, and lU is a 
normalization constant to make fz^, (z) dz = 1. Besides, we assume Rayleigh fading in 
this subsection since Rayleigh fading is the worst case for the proposed analysis as addressed 
in Subsection IIV-AII 

In the following, we investigate the considered network with the proposed microscopic 
analysis of the UE interference distribution. Eirst, for each b E {2,..., 84}, we invoke 
([9|), mH) and ([TSl) to check the maximum error e among be {2,..., 84} using ([20l) . If the 
maximum value of e is reasonably small, e.g., less than 0.01, then we can approximate 
in ([5]) as a power lognormal RV = 10^*^, where the PDE and the CDE of Q are given 
by ([2T[) with the parameters A, pq and Gq obtained from, e.g.. Appendix B. 
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The maximum values of the 83 i?b-speeifie e’s for various r values are presented in Table 0 
From Table 0 we ean observe that, for all the investigated values of r, the maximum values 
of £ are below 0.01. Thus, each Ib should be well approximated by a Gaussian RV Qb- Due to 
space limitation, we omit the detailed numerical investigation on the Gaussian approximation 
for each Ib, which is very similar to the discussion in Subsection IIV-A2I After obtaining the 
approximation for each Ib, we approximate in (|5]) as a power lognormal RV 7™"^ = 
using (l22l) . The numerical results of A, hq and ctq are provided in Table 0 for reference. 


Table V 

Approximation errors of the proposed analysis {B = 84, non-uniformly distributed UEs with the 
DISTRIBUTION FUNCTION OF /z^ (z) = ^,Z £ Rb, RAYLEIGH FADING). 


r 

Max e 

Max El + E2 

Max El 

Max s[ + e '2 

Max e'l 

A 

MQ 


0.01 

5.1 X 10-3 

2.7 X 10-4 

4 X 10-3 

4.8 X 10-3 

4 X 10-3 

48.9 

-99.7 

116.2 

0.02 

5.1 X 10-3 

4.2 X 10-4 

4 X 10-3 

4.7 X 10-3 

4 X 10-3 

48.0 

-101.6 

116.5 

0.04 

5.9 X 10-3 

1.2 X 10-3 

4 X 10-3 

4.7 X 10-3 

4 X 10-3 

47.5 

-103.1 

117.2 



Fig. 6. The simulation and the approximation of the CDF of 1"’"^ in dBm {B = 84, non-uniformly distributed UEs with 
the distribution function of /z^ {z) = ^, z £ Rb, Rayleigh fading). 

To further verify the accuracy of our analytical results on the UL interference distribution, 
in Fig. we plot the simulation results of the CDF of in dBm and the approximate 
analytical results according to (l22l) . As can be seen from Fig. |6l the resulting power lognormal 
approximation of is tight. However, note that the approximation shown in Fig. is 
not as perfect as that exhibited in Fig. [3l According to the discussion in Section |nll the 
approximation errors associated with the first and the second steps of approximation are 
captured in e, which is very small as can be confirmed from Table 0 The noticeable small 
approximation errors in Fig. are caused by the inaccuracy of approximating the sum of 
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multiple lognormal RVs as a single power lognormal RV. Note that in our previous work ifTOll . 
we use a lognormal RV to approximate the sum of multiple lognormal RVs, which leads to 
larger errors compared with the results shown in Fig. 0 Finding an even better distribution 
to approximate the CDF of the sum of multiple lognormal RVs will be our future work. 

C. Discussion on the Complexity of the Proposed Microscopic Analysis 

The computational complexity of the proposed approach is mainly attributable to the 
numerical integration required to obtain the values of e for each small cell. In contrast, 
the simulation approach, e.g., lUl, BUl, as well as in this paper, involves a tremendously high 
complexity. Specifically, in our simulations, in order to go through the randomness of all the 
RVs discussed in Section [III more than one billion of realizations of Ib have been conducted 
for the 83 interfering cells depicted in Fig. [51 This shows that the proposed microscopic 
analysis of network performance is computationally efficient, which makes it a convenient 
tool to study future 5G systems with general and dense small cell deployments. Furthermore, 
our analytical studies yield better insight into the performance of the system compared with 
simulations. 


V. Conclusion 

The lognormal approximation of the UL inter-cell interference in FDMA SCNs is important 
because it allows tractable network performance analysis. Compared with the existing works, 
in this work we have analytically derived an upper bound on the error of such approximation, 
measured by the KS distance between the actual CDF and the approximate CDF. 

Our results are very general in the sense that we do not pose any requirement on (i) the 
shape and/or size of cell coverage areas, (ii) the uniformity of UE distribution, and (iii) 
the type of multi-path fading. Based on our results, we have proposed a new approach to 
directly and analytically investigate a complex network with practical deployment of multiple 
BSs placed at irregular locations, using the approximation of the aggregate UL interference 
by a power lognormal distribution. From our theoretical analysis and simulation results, we 
can see that the proposed approach possesses the following merits: 
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1) It quantifies the approximation error measured by an upper-bound KS distance using 
a closed-form function. And the tightness of the approximation is validated by the 
numerical results. 

2) It tolerates more practical assumptions than the existing works, e.g., irregular hot-spots, 
overlapped cells, etc. And it can cope with a large number of small cells with a low 
computational complexity of analysis, thus making it a convenient tool to study future 
5G systems with general and dense small cell deployments. 

As future work, we will further investigate the impact of the correlated shadow fading, the 
three-dimensional (3D) antenna pattern and the multi-antenna transmission on the proposed 
approximation of the UL interference distribution. 

Appendix A: Proof of Theorem [2] 

For clarity, we first summarize our approach to prove Theorem [2] as follows. Our idea is 
to perform a Fourier series expansion for both the CDF of G and that of the hypothetically 
approximate Gaussian RV. The distance between those two CDFs will be quantified by the 
upper-bound KS distance derived in closed-form expressions. 

First, according to the definition of G, the CDF of G can be formally represented by 


Fg (^) = Pr G <g 


= Pr 


= Pr 


L + S <g 
S<g-L 


r’ + OO 


$ 


9 - I 
(^s 


h (0 di, 


(23) 


where $ (■) is the CDF of the standard normal distribution. According to [fT^ . $ can 

be further written as = 1 ~ ^crfc where erfc (x) is the complementary error 

function defined as erfc (x) = exp (—t^) dt in ifTSll . 

Second, due to the independence of L and S, the variance of the approximate Gaussian 
RV of (L + S) should he = 0 and cr£ + cr|, respectively. As a result, the CDF of 


(24) 


such approximate Gaussian RV can be expressed as 


$ 


9 




= 1-erfc 

2 


^2 (al + al) I 
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In the following, we will derive the elosed-form expressions of the upper-bound KS distanee 
between the two CDFs presented in (l23l) and (l24l) . respeetively. 


(25) 


According to [l22ll . erfc (x) can be expanded to a Fourier series as 

A —1 / 9 9\ 

^ / \ 4 exp [—n oj ) , . ^ , 

erfc (x) = 1-> -sin iznuxj -f d (x), 

IT n 


n=l,n odd 

where oj is the fundamental frequency of the Fourier series, N = 2p—1 is the series truncation 
point, and 5 (x) is the residual error of the p-truncated Fourier series. 

Based on (|25]) . ^ in (1231) can be expanded as 


$ 


9 - I 
0-5 


2p-l 


^2 + x ^ 

n=l,n odd 


exp {—rroj^) 


n 


1 2 exp{-n^u^) 

9 TT ^ 


2 n 

n=l,n odd 


n 


- Is ( 

V V 2 asJ 2 \y/ 2 asJ 

imag |exp(j^/2»a»^)} - Is , (26) 


where imag {■} extracts the imaginary part of a complex value. Plugging (|2^ into (|2^ . yields 

OnMnX PYn t/■+°° / \/Oni,)i\ 


2. I f .y/2noj(^exp (—n^u'^) f~^°° ( .^/2nul^ 

> exp j- - 

f ^ \ fTn I n 

lj 2 =l,n odd 


-imag 

TT 1 ^ " V" 


/ exp -j- 

oo \ 


where the integral exp j fi{l) dl is actually the characteristic function ITT^ of 

fl (1) evaluated at the point t = — Such characteristic function, denoted by ipi (t), is 

^-1-00 

(t) = / exp (jtl) fl (/) dl. 


(28) 


Hence, from (1271) and (1281) . (1231) can be re-formulated as 

(ixl'] 

\V2asJ 

j 


1 1 /■+°° 


f fl (/) dl 


2p-l 


2. I /. 2/2nojg \ 1 / 2 9\ / V^nu 

-imag > exp j- - exp [-n u))^i 

Ln=l,nodd \ a / \ a 


r *-\-00 


( ^T^jh (0 dl + -imag \ V exp (j 


f . '/2nujg \ 

^ A 


where Vr, is defined as 


A 1 / 2 2 \ ( V^nu 

Vn = - exp i-n OJ )ipi\ - 

n ' \ Os 


(29) 


(30) 
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And the summation in the imag {■} operator of (l29l) ean be deemed as the weighted sum of 
the unit-amplitude complex values jexp n G {1,3,...,2p — 1}. 

Regarding Vn, by means of the series representation of exp (•), we can re-write it as 

1 / 2 2\ I 

Vn = — exp (^—n u ) (pi 


n 


(ys 


= - exp 

n ^ ' 


= - exp (-n^uj^) 
n ^ ' 


^+oo 


exp 


^+oo 


1-j 


<ys J 

.\/2nul 1 f y/2nul\ 


. 1 f y/2nul\ 1 f y/2nojl\ 


as 2! 

4 


\ (ys J 
h (0 dl 


1-0 


= — exp (— 
n ^ ' 

(“) 1 ( 2 2\ ( 

= — exp ( —n u 1 exp 

n ^ ' \ 


2n^u‘^ar . .pf'^ (\f2‘ 


2a2 6 


.(4) 


nu \ Pl' I \/2 nu 


(ys 


24 \ as 


n^u'^al 


ox 


+ f^r. 


= — exp —n u 
n 


2, ,2 / '^L + 


Ot 


+ ly-v 


(31) 


where v„ is defined as 


2 21 (yf + (y^ 

Vn = - exp I -n u ' 


A 1 


n 


ot 


(32) 


Besides, in (l3TI) . measures the difference between Vn and defined as Hn = Vn — Vn- 

The motivation of introducing Vn into (l3TI) is because (1 — ^ ) can be approximated as 

( 2 2 2 \ 2 2 2 _ 

) when is small in the step (a) of (ISTl) . Note that no approximation is 

'^s J 

assumed in (l3TI) because fully captures the difference between Vn and Vn- 


Plugging (I3TI) and (l3^ into (l29l) . we can get 

. ' 


2p-l 


1 ^=— j fi (/) dl + -imag ^ «n exp i 

^ fn=l,nodd \ 


2 1 


TT 


n=l,n odd 


exp I + Og ^ ^ 
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Performing a variable change of ce = ce 


i in (1^ . we can obtain 


1 1 /■+°° 




exp 


i-\/2nojg 

^/ol + Vl 


9 2p-l / 

i 5: iexp(-„V)=i„ - 

n=l,n odd \ 


\/2nug 


y^i + ^s 

On the other hand, from (|25]) . (I2^ can also be expanded to a Fourier series using Cj as 


(34) 
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^1-1^ 


a/2 ((t£ + 


CTc 


2 1 


i E iexp(-nV)si 


sm 


TT ' n 

n=l,n odd 


\/2nug 


(35) 


where 5 (•) is the residual error function incurred from the p-truncated Fourier series expansion 
and it has been defined in (|25]) . 


From (l34l) and (1351) . we can bound the distance between Fg (g) and <h ( —r=—; 
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(36) 


n=l,n odd 

The first two terms in the right-hand side of (l36l) are caused by the residual errors from 

the p-truncated Fourier series expansion. From ll23]l . 6 (x) in (l25l) can be strictly bounded by 

/ (a:)| < ^ erfc ((2p + 1) a;) + erfc (—— bl) . (37) 

V TTo; V 20; / 


Note that such bound is dependent on the value of x, which makes it difficult to handle the 
first term in the right-hand side of (l3^ . because the integral in that term is computed from 
—cxo to -l-cxo. Thus, in the following, we will discuss the upper-bound of / (a;)| respectively 
for the case of |x| < A: and for the case of |x| > k, where A: is a positive scalar. 
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Considering the monotonic increase of erfc — |a;|) with respect to |a;|, we can further 
bound (x)| for a range of interest of |x| < A: as 

|(5 (x)| < ^^erfc ((2p + 1) te) + erfc - A:) = (w,p, A:), (|x| < A:). (38) 

■yTTCJ \2UJ / 


For the range of |x| > k, since erfc (x) < 2, |(5 (x)| can be bounded by 

|() (x)| < ^^erfc ((2p + 1) w) + 2 = (a;,p), (|x| > k) 


(39) 


In order to derive an upper bound that is independent of g for the right-hand side of (l3^ . 
in the following we consider two cases for g, i.e., < ki^/+ cr| and \g\ > A;iy^o^^T”^, 

where ki is an arbitrary positive scalar. 

When \g\ < ki^^al + cr|, we can bound the first term in the right-hand side of (l3^ by 
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where k 2 is another arbitrary positive scalar introduced to facilitate the bounding of the 
to -foo with ] 
f-L (0 dl = Pr 


integral from — oo to -foo with regard to 1. The last step of (|4^ is valid because (i) 
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> k2\lal + al 


< Pr 
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< (41) 
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which comes from Chebyshev’s inequality Il2^ . and (ii) f hil) dl<Jfi{l) dl = l. 

\l\<k2y/al+al ^6® 

Besides, the second term in the right-hand side of (l36l) can be bounded by 
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V2f 2 

Plugging (l40l) and (1^ into (l3^ . and considering the definition of £i and £2 in (|9l) and (fTTI) . 
respectively, we can get 

(a') - I ; / . j <£i +£2, ( \g\ < ki\l(j'i + (y's ] ■ (43) 


— ^ <£i + £ 2 , [\g\<ki\lal + a\ 

V^L + 

When \g\ > ki^af + cr|, we invoke Chebyshev’s inequality [|24]| to obtain 
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Pr 
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> kiGal + al 
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Therefore, we ean bound 
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where (l45l) and 0 < 1 — <h , —; 

into the seeond last step of (l46l) . 


< |erfe (A:i) for | 5 '| > ki^yaj^+~^, have been plugged 


Combining (1431) and (l46l) . and eonsidering the definition of £3 in (fT^ . we ean obtain ([ 8 ]), 
whieh eoneludes our proof. 

Appendix B: An Example to Obtain A, and gq 


Aeeording to [EU, with regard to the RV Q, we have 

q^+oo oq UQ 

d^-^[FQ{q)) ^ ^ 

dq o-Q 


(47) 
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q^—oo 


Besides, aeeording to ifTTll and [|20ll . we ean get 
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(48) 


where ax is obtained by solving the following equation set IfTTll 

(si) = n (■si) = Cl 

b=2 

(S2) = n ^Qb (^2) = G2 

b=2 


(49) 


where (s) is the approximate moment generating funetion (MGF) evaluated at s for a 


lognormal RV defined as 10 10 ^. Sueh approximate MGF is formulated as 
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is) = 


Mo / 

E 'U’m I 

T^exp -.exp 

m=l ^ \ 


V2 


“ 1 “ /^X 


(50) 


where C = Mq is the order of the Gauss-Hermite numerieal integration, the weights 
Wm and abseissas am for Mq up to 20 are tabulated in Table 25.10 in If251 . Usually, Mq is 
set to be larger than 8 IfTTlI . Similarly, in (l49l) . (s) is eomputed by replaeing fix and 

ax respeetively with fiq^^ and aq^^ in (l50l) . In (l49l) . si and .2 are two design parameters 
for generating two equations that ean determine the appropriate values of fix and a\. For 
example, we ean choose si = 0.001, .2 = 0.005 and Mq = 12 as recommended in [[T71l . The 
solution of (l49l) can be readily found by standard mathematical software programs such as 
MATLAB. Besides, using fix obtained by solving (|4^ . we can match the mean of Q with 
fix to construct a third equation to determine the three parameters, i.e.. A, fiq and aq, in the 
power lognormal distribution of Q. 

To sum up, based on (l47l) . (|4^ and matching the mean of Q with fix, we have the following 
equation set to determine the values of A, fiq and aq. 


1 _ 1 

(a) 


O-Q (TX 


vA _ rr~‘^ 

UQ - \J Mb=2 

(b) ■ 

(51) 

<lfQ (^) 

(c) 



Note that equations (a) and (b) of (I5TI) are easy to solve and they can deliver A and aq, 
while equations (c) of (l5TI) can be efficiently solved by the standard bisection process |[2^ 
to determine fiq. 
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